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The current status of the theoretical understanding of the transverse-momentum dependent 
. parton densities (TMDs) is discussed. Special attention is payed to the difference between 

the operator definitions of TMDs proposed so far, the treatment of specific divergences, 
^ ^. the geometry of the gauge links, and the role of the soft factors. 

(D ■ 

(— I i Short-distance factorization is the basic concept for the application of QCD at high energies 

[U [21 13] . The predictive power of QCD in this region has been successfully demonstrated in the 
case of fully inclusive processes, where the coUincar (Feynman) parton distribution functions 
I (PDFs) arc applicable as a nonperturbativc part of the QCD factorization theorems. Appli- 

^SJ ■ cation of the QCD factorization approach to the description of the semi-inclusive hadronic 

0^ ' reactions, such as the scmi-inclusivc DIS, or the Drell-Yan (DY) process, demands an expan- 

. sion of the concept of parton distribution functions beyond the coUincar approximation. It is 

useful to introduce the transverse-momentum-dependent (TMD) parton distribution (or frag- 
04 . mentation) functions (in what follows I abbreviate them by TMD = transverse-momentum 

' "densities") are essentially non-pcrturbative objects, which accumulate information about the 

intrinsic 3-dimensional motion of partons (quarks and/or gluons) in a hadron. In Refs. [HEIEI, 
the TMDs have been proposed as a generalization of the collinear PDFs: 

; -P[collinoar](a;,M) ^ -^[tmd] (a^, kj_ , /i, C) • (1) 

H 

■ It was pointed out that, in contrast to the collinear case, TMD depends on "how fast the hadron 

in moving" , that is, on hadron's rapidity. Formally, this dependence enters the Eq. ([1]) via the 
additional variable C = (2P • v)/\v'^\ , v'^ ^ 0, where P is the incoming nucleon momentum 
and 71 is non-light-like axial gauge fixing vector. The appearance of the additional rapidity 
dependence in the TMDs requires the development of an appropriate rcsummation scheme and 
the derivation of suitable evolution equations. Another important consequence of the rapidity 
dependence that has been made by Soper [4]: the very definition of a parton is what changes 
as rapidity varies. In the collinear case, the only extra variable is the rcnormalization scale 
/i, and the gauge-invariant operator definition of integrated PDFs allows one to relate their 
moments to the matrix elements of the twist-two operators arising in the operator product 
expansion on the light-cone [7] , thus making the PDFs consistent from the point of view of the 
local quantum field theory and allowing the probabilistic interpretation in terms of the parton 
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number operators. The renormalization properties of the colhnear PDFs are governed by the 
DGLAP equation, estabhshing the logarithmic dependence of the structure functions on the 
hard scale . 

Need of a consistent theory of the TMDs as quantum objects initiated the quest for the 
deeper understanding of such fundamental issues as QCD factorization, universality of PDFs, 
their renormalizability and complete gauge invariancc. The next important step in the theory 
of TMDs had been done in the works [71 |H]: it was observed that trying to write down the 
operator definition of the (quark) TMD in the form similar to the colhnear PDF (the axial 
gauge {A - v) ~ , ^ is used to remove the longitudinal gauge links) 

^i/h(x,k^,M,C) = \j ttX^^ '"''^^'^'''^^^ {h\m-.i±)l^UQ-.^±)\h) , (2) 

where the extra parameter is defined via the gauge-fixing vector C = (2P • w)^/|t;^|, one is able 
to separate out the emergent rapidity divergences in the form of the (powers of) logarithms 
In ^ and perform their resummation by means of suitable evolution equation [7] . Although the 
formal integration of definition ^ over yields the colhnear PDF 

j (fikx_ J^i/h(a;, k±, /X, C) -F[coninoar](a;, m) , (3) 

such a transition can only be justified in the tree approximation, since the presence of the 
rapidity divergences breaks down this direct procedure. Nevertheless, important progress has 
been achieved in the semi-inclusive factorization (in particular, for the Drell-Yan process) in 
the works [3 [5] and [S] • 

Since the beginning of the nineties, a wide range of the phenomenological application of the 
TMDs has been revealed in the polarized hadronic processes. In particular, it was found that 
TMDs play a crucial role in the understanding of the single-spin asymmetries (see, e.g., Ref. 
[ini HH [m [El [HI m] and Refs. therein). It has been realized that a much deeper knowledge of 
the properties of the polarized and unpolarized TMDs is essential not only because of the great 
theoretical interests in them, but from pragmatic needs: one has to be able to use the correct 
operator definitions of TMDs, to know the complete set of their evolution equations, to have 
justified factorization theorems, to keep under control the universality of the TMDs, etc. Since 
then, important steps have been done in the study of the initial and final state interactions 
in the operator formulation of TMDs and their relations with the structure of Wilson lines 
[TBI [T71 [T5] ; TMD factorization in the covariant gauges [THj ; treatment of the extra singularities 
within the subtraction scheme j20[ I21j ; problems of the universality and gauge invariance of 
TMDs and structure of the gauge links [22l [23] ; the problem of (breakdown of) factorization 
in certain semi- inclusive processes [Ml [2S1 [26] ; problems related to the emergent self-energy 
singularities of the Wilson lines entering the definition of TMDs [27]; the issue of matching 
different methods of taking into account the kj^ -dependence at different momentum scales [28) : 
various aspects of the evolution equations for the TMDs [351 [30]; calculations of the TMD 
within the quark models [31], development of the lattice simulations methods for the TMDs 
[37] . development a generalized approach to the operator definition of the TMDs taking into 
account their renormalization-group properties in the light-cone gauge [34 1 135 1 136] . stc. Further 
important references can be found, e.g., in Refs. [32 1 133 j . 

The main problem I intend to discuss in the present work is the treatment of pathological 
singularities beyond the tree-approximation in the TMDs. To be precise, at the one-loop level 
the following sorts of singular terms arise: 
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• Ultraviolet poles - in the dimensional regularization: they have to be removed by the 
standard renormahzation procedure. 

• Rapidity divergences: they depend on the additional rapidity parameter C [3 13 ; but do 
not break the renormalizability of the TMDs, and can be safely resummed by means of 
the CoUins-Soper equation. 

• Pathological overlapping divergences: they contain the UV and rapidity poles simulta- 
neously ^ - In 77 . They break down standard renormalizability of TMDs, calling for a 
generalized renormahzation procedure in order to enable the construction of consistent 
operator definition of the TMDs. 

To overpass the problem of emergent extra divergences, different approaches to the operator 
definition of TMDs have been proposed. In what follows, my notations are: A or C denote 
the axial or covariant gauge, respectively, and the subscripts v or n mark the vectors defining 
the non-light-like or light-like directions of the longitudinal gauge links (in covariant gauges) or 
the gauge- fixing vector (in axial gauges), respectively. I will distinguish between the following 
definitions of TMDs: 

• Av-TMD: one applies the axial non- light-like gauge, the longitudinal (along the vector 
v) gauge links vanish by virtue that {v ■ A) = 0; the rapidity cutoff is defined as C = 
{2P-v)^/yi Refs. [ilH]. 

• Cv-TMD: in the covariant gauge, the longitudinal gauge links along the non-light-like 
vector survive, the transverse gauge links cancel; the rapidity cutoff ^ = (2P • v)'^/\v'^\ 
characterizes the deviation of the longitudinal gauge links from the light-like direction; 
the soft factor contains the non-light-like gauge links, Refs. jl9j . 

• An-TMD: one uses the light-like axial gauge {n ■ A) ~ , ~ Q, the light-like longi- 
tudinal gauge links vanish, the transverse gauge links at the light-cone infinity survive; 
the regularization parameter is defined as 77LC = (^ ' '^)/'7i where 7] originates from the 
regularization of the g^-pole in the gluon propagator, e.g., with the PV prescription 

11/1 1 

[q+]rj 2 \q+ + iij q+ - ir] 

the soft factor contains the light-like longitudinal and the transverse gauge links, Refs. 

m- 

• Cn-TMD: covariant gauge with light-like gauge links, no transverse gauge links; soft factor 
contains both light-like and non-light-like longitudinal, as well as the transverse gauge 
links; regularization parameter C, = (2P ■ w)^/|w^|, where v is the non-light-like direction 
of gauge links in the soft factor, Refs. [20 l [2T 1 [22 | |32] . 

• L-TMD: this definition is proposed for the lattice simulations; the direct connector is used 
as a gauge link, no regularization parameters, no light-like gauge links |37j . 

It is worth noting that there is no a priori direct relationship between these definitions, they 
define, in principle, different objects. 

Let us discuss the properties of the emergent singularities within the definitions listed above. 
Within the "off-light-cone" frameworks Av-TMD and Cy-TMD, a rapidity divergence arises in 
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the form of the logarithmic terms In^ and In^ ^. In the "hght-cone" definitions An-TMD and 
Cn-TMD, the situation differs from the previous case. It was shown in [34] that in the An- 
TMD. the anomalous divergent term containing overlapping (UV (g) rapidity singularity) stems 
from the virtual-gluon contribution (we consider the "distribution of a quark in a quark with 
momentum p" and the TMD is normalized as 7r'^"''°''(a;, kj^) = S{1 — (k^)) 



S[A„] 



-Cf r(e) 



^ [^\v 



(5) 



where the x- integral is regularized according to the rule In the covariant gauge, the defini- 
tion C„-TMD yields [S] 



-CFr(e) 



4n- 



5{l - x)5^^Hk^) 



dx 



(I-x) 



l+e 



(6) 



One easily observes that without the Ty-regularization of the last integral and after a trivial 
change of variables, Eq. ([5]) is reduced to Eq. I][Cn](e) ~ S[An](e,7; = 0) . In principle, one 
can use dimensional regularization to take care of the overlapping singularities as well. I will not 
discuss here this approach. In more detail, the applicability of dimensional regularization for a 
consistent treatment of the divergences arising in the path-dependent gauge invariant 2-quark 
correlation functions had been studied in Ref. }38j . 

Within the An and Cn , the cancellation of the overlapping divergences is achieved by means 
of the subtraction of the special soft factors. This procedure implies the generalized renor- 
malization of the TMDs and provides us with a completely gauge invariant object, free of the 
overlapping divergences [341 135j . However, another problem arises in the one-loop corrections 
to the soft factor itself. In the light-cone gauge, one obtains the term 



Ssoft[A„] =zgV'CF2p+ J^"^ 



1 



(27r)" q^{q- ■ p+ - iO)[q^ 



(7) 



where a new singularity emerges, that can not be removed neither by dimensional regularization, 
nor by the 77-cutoff. One has 



■^soft 



[An] 



-Cf 



47r/i 



21 « 



A2 



r(e) 



dx- 



(8) 



where A is the IR regulator. Taking into account that the extra singularity is cusp (or rapidity) 
-independent [34] . we conclude that it represents the self-energy of the Wilson line evaluated 
along a "straightened" path, i.e., with the angle becoming very small: P+ — rj. Subtraction of 
this self-energy part is presented graphically in Fig. 1. Note that there is no need to introduce 
additional parameters in this subtraction and this procedure doesn't affect the rapidity evolution 
equations. Moreover, it has a clear physical interpretation: only an irrelevant contribution due 
to the self-energy of the light-like gauge links is removed, which is only an artifact of the 
unobservable background. 

Therefore, the following generalized definition of the An-TMD is proposed [34l ESI [36] 



(a;,kj_;pi,r/) 



d^d^U 



_l; 00 



i/h v^,"^,^,.,; - 2 J 2^(2^)2 
;[oo~,C_L;oo",oo_L]'l'7+[oo",oOi;oo",0_L][oo",0j_;0~,0^]'i/;i(0",0_L)|/i) i?"^ 



(9) 
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Figure 1: Subtraction of the Wilson- line self-energy contribution in the soft factor. 
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where the gauge invariance is ensured by the inserted gauge links [y, x\ =1^ exp \_—ig /J dz^A^{z)\^ 
with A = t"'A'^, and the contours for the soft factor are presented in Fig. 1. 

It is worth mentioning that within the An-TMD, the soft factor can be interpreted in terms 
of the "intrinsic" phase resembling the Coulomb phase found by Jakob and Stefanis in Ref. j39] 
in QED. Its origin was ascribed in [39j to the long-range interaction of the charged particle with 
its oppositely charged counterpart that was removed "behind the moon" after their primordial 
separation. Such a reminiscence of the quantum entanglement has been used recently |40j in 
our analysis of the long-rapidity-range two-particle correlations produced at the LHC (reported 
by the CMS collaboration [41], see also these Proceedings). Therefore, within the An-TMD 
approach, one may understand the appearance of the soft factors in the TMD factorization not 
only in technical terms, but also from the basic principles of a quantum theory. 
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